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We propose a new formulation for atomic side mode dynamics from super-radiant light scattering 
of trapped atoms. A detailed analysis of the recently observed super-radiant light scattering from 
trapped bose gases [S. Inouye et al., Science 285, 571 (1999)] is presented. We find that scattered 
light intensity can exhibit both oscillatory and exponential growth behaviors depending on densities, 
pump pulse characteristics, temperatures, and geometric shapes of trapped gas samples. The total 
photon scattering rate as well as the accompanied matter wave amplification depends explicitly on 
. . . atom number fiuctuations in the condensate. Our formulation allows for natural and transparent 

• interpretations of subtle features in the MIT data, and provides numerical simulations in good 

I agreement with all aspects of the experimental observations. 
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I. INTRODUCTION 



The successful discovery of Bose-Einstein condensation 
(BEC) in dilute trapped atoms provided significant 
momentum for research into quantum degenerate gases. 
In analogy with laser theory, condensation results in a co- 
' herent matter wave field, which has since been identified 
and several important optical analogous effects in- 
eluding four- wave mixing |^ , superradiance , coherent 
[ matter wave growth |^ were demonstrated. Theoretical 
studies of these phenomenon in degenerate BEC systems 
[^,0 pointed out the important role of correlations and 
competitions among matter wave side modes, i.e. multi- 
Ch ' mode nature of even a single component condensate due 
. to center of mass (CM) motional effects. For a harmoni- 
^ ■ cally trapped condensate, these multi-modes can be con- 
Q , veniently expressed in terms of quantized motional states 
O ' with equal energy spacing. Theoretical investigations of 
^ \ light scattering from such trapped atoms are complicated 
since both the elastic and inelastic spectra can include 
contributions from many different motional states. 
, In this article, by proposing a new identification of 
■ - ' trapped atomic side modes for light scattering from a 
plane wave excitation, we attempt for a detailed interpre- 
tation of the recently observed off-resonant super-radiant 
light scattering |^ . This constitutes an example of using 
light scattering as a spectroscopic tool to probe proper- 
ties of trapped degenerate quantum gas. Our investiga- 
tions show that quantum statistics of the condensate can 
have a drastic effect on the properties of scattered pho- 
tons |^,||. Our formalism takes advantage of the recent 
success with atomic multi-modes pO| to provide a clean 
interpretation of all aspects of the experimental observa- 
tions jij. Similar approaches can also be used to clarify 
physical pictures of the more recent BEC Bragg spec- 
troscopy experiments jl^ . 
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II. FORMULATION 



Our model describes light scattering of trapped atoms 
from excitation due to an intense far off-resonant plane 
wave pump field. The proto-type system is illustrated 
in Fig. |l| as arranged in the MIT experiment [Q. The 
atomic sample is assumed dilute and atoms are of the 
alkali type with a single valance electron. Two electronic 
states [g (e) for ground (excited)] are connected by a real 
electronic dipole moment d. In units of fi. = 1 and in 
length gauge, our model Hamiltonian takes the form 

a=g,e-' 

+ il7o(t) / dr-^l{r)e'^'>-^-^g{r) + h.c. 



+ dr dkg{k)-i>l{r)e'^%-fg{r) + h.c., (1) 



under dipole and rotating wave approximation. Ha^ — 
—'V'^ /2M+Va{R)+Ecr is the atomic Hamiltonian consists 
of CM kinetic energy (atomic mass M), trapping poten- 
tial Va, and electronic excitation energies in the rotating 
frame Eg = 0, E,, = A = uja — ojq. A is large since the 
pump field at frequency ujq is far detuned from the atomic 
transition frequency uja- The dipole interaction between 
an atom and the pump field is described by a time de- 
pendent Rabi frequency ^o{t) — QoT{jot), with Qq the 
peak value and 70 the temporal width of the envelope 
function T{jot). Both fio and T{'yot) are chosen to be 
real. The polarization and wave vector of the pump pulse 
are denoted by eo and ko, respectively. The second term 
of Eq. (|l|) describes free electromagnetic fields needed to 
consider inelastic photon scattering. Their polarization 
index is suppressed in the wave vector k. The degenerate 
atomic fields are described by annihilation (creation) op- 
erators ^(r) [^'^(r)] satisfying appropriate commutators 
for bosonic or classical Maxwell Boltzman statistics. The 
operator 6^ (6t) annihilates (creates) a photon with wave 

vector fc, polarization e^, and energy LUk — cfc — (again 
in the frame rotating with frequency luq). Within the two 
state approximation, dipole coupling between scattered 
field and an atom is given by g[k) — —iy^ ckd"^ /An'^e^ ■ e^:, 
a slowly varying function of k. When the detuning A is 
much larger than any other frequency scale in the sys- 
tem, the excited state can be eliminated to obtain an 
effective field theory within the ground atomic motional 
state manifold ||], 

H^J dr^lif) \^HAg + *,(f) + I dku,b\bj: 

+ [ df [ dk g{k,t)^l{r)e'^''-^"'^-%'^g{r) + h.c. . (2) 



The AC stark shift term from all field modes except the 
pump will be neglected in this study since we are in- 
terested in regimes when the scattered field intensity re- 
mains small. The effective coupling constant g{k,t) — 
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f^o^(7o05(^)/2A is now time dependent, and it de- 
scribes the scattering event in the ground motional state 
manifold, i.e. the combined effect of absorbing a pump 
photon to the excited state manifold followed by a spon- 
taneous emission back to the ground state again. In most 
cases g{k, t) is a slowly varying function of time. To sim- 
plify our discussion we consider a temporal square shaped 
pulse and thus drop the constant term |rio(^)P/2A. In 
contrast to strongly correlated many body systems the 
trapped atoms are assumed non-interacting, allowing 
for detailed investigations of their interaction with light 
fields. The ground state atomic field operator can be ex- 
panded in terms of trapped single atom wave function (/>„ 
according to g{r) — ^„ c„0„(r). For a harmonic trap- 
ping potential, 0n(r) is simply the number basis states in 
position representation with n = {rix, ny,nz) a triplet in- 
dex. Atomic annihilation and creation operators c„ and 
cl obey bosonic algebra [c„, cJ,J = dnrn and [c„, c,„] = 0. 
We then reduce the Hamiltonian to 

H = Ho + Hr + Has + Hs , 
Hq = '^EnC^Cn / dkuJkbibj:, 

n 

Hr^^ dkg*{k,t)rin,n{k- ko)clbicn + h.c., 

Has ^ dkg*{k,t)i]n,m{k- ko)clbic,n + h.c., 

ffs = ^ ^ dkg*{k,t)r]n^,n{k- ko)clbic,n + h.c. . (3) 

Figure ^ is a pictorial display of characteristic absorption 
and emission cycles for Rayleigh, Raman Stokes and anti- 
Stokes processes. The atomic energy is = -Eq + w„, 
consisting of contributions from electronic ground state 
energy Eq — and CM motion energy cun — i-^t ■ n 
with frequency cJt = {^xi^yT'-^z) for a three dimensional 
trap. The factor rin,m{K) = {n\exp {—iK ■ f^\m) repre- 
sents CM motional state dipole transition moment, and 
is analogous to the Franck-Condon factor in a di-atomic 
molecular transition. It is simply the matrix element 
of displacement operator D{r) = exp (— • r) in the 
number basis and depends on the total recoil momentum 
K = ko — k from the scattering cycle involving absorp- 
tion of a pump photon followed by an emission. Within 
the ground motional states, it acts like a diffraction ma- 
trix since it shifts atomic fields around in momentum 
space. To examine various competing dynamical pro- 
cesses in light scattering described by Eq. (^) we sepa- 
rate the coupling term depending on the energies of the 
two coupled (initial and final) motional states. This leads 
to three types of scattering (as in Fig. |2|): 1) the elas- 
tic Rayleigh scattering described by Hr corresponds to 
events within the same atomic motional states; 2) the 
Stokes Hs and 3) the anti-Stokes Has Raman scattering 
into higher (lower) energy motional states. We empha- 
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sise that Stokes and anti-Stokes terms here corresponds 
to the same final electronic state but with increased or 
decreased energy final motional states. We may thus also 
call them as "inelastic" Rayleigh scattering. 

Before a detailed discussion of the three type scattering 
events, it is possible to get a crude picture of how each 
individual interaction term contributes to the dynamics 
at low temperatures and within a short time scale. When 
the gas sample is at sufficient low temperatures only low 
lying atomic motional states are densely occupied and 
therefore their atomic fields can be approximated as clas- 
sical variables, while upper motional states are sparsely 
populated and need to be treated quantum mechanically. 
Within a short duration, an approximation involving un- 
depleted populations in lower motional states, similar to 
the parametric pump approximation in nonlinear optical 
multi-mode coupling, can be made. We can then consider 
a single motional state (m) and a single resonant scatter- 
ing field mode with k — — ko), and assume all A^o 
atoms were initially condensed in state n = (0, 0, 0). We 
then approximate cq ^ Cq ^ VNq, and Hamiltonian (j^) 
is simplified to 

H^Ho + Hr + Has + Hs. 

Hr = 9*{ki,t)r]a^o{ki - k())Nabi^ + 0{cl^c„,bf:J + h.c, 
Has = g*{ki,t)i]o,m{ki - ko)\/l%b\ Cm + h.c, 
Hs = <?*(fci,t)?7™,o(fci - fco)V^cJ„&^^ + h.c. . (4) 

With this simplified model, Hr becomes a displacement 
operation for scattered field mode (fci); it describes the 
generation of coherent photon fields in this mode. For the 
general case of a thermally populated motional state dis- 
tribution, the total coherent photon fields are distributed 
accordingly. Terms of order 0{cly^Cmb^^) are neglected 
since the second order quantum processes described by 
Hs and Has are more important within the short time 
period discussed here. We see that Has resembles an 
atom-polariton Hamiltonian, describing the generation of 
atom-photon bound states; while Hs takes the form of 
a non-degenerate parametric amplifier Hamiltonian, de- 
scribing processes of gain, squeezing, and atom-photon 
entanglement. It should be noted, however, even when 
both terms Hg and Has are in existence, they can still 
be grouped as a more general polariton Hamiltonian. 
Thus, internal fields, i.e. scattered photons inside the 
atomic sample can always be viewed as being part of an 
atom-polariton system Overall photon scattering 

and emissions from such a system is complicated as in- 
ternal fields couple to external fields to cause radiative 
decays of atom-photon polaritons. Their periodic energy 
exchange implied by Has could exhibit Rabi type oscil- 
lations in the radiated fields. If on the other hand when 
the motional state m is initially unpopulated, oscillatory 
behaviors will not appear since Has can not contribute 
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during the early stages of dynamic interaction. 

In the foUowing two sections, we wiU focus on situ- 
ations when Hs is the dominant interaction term. We 
start by outhning the required system conditions to 
achieve this. First we observe that Hn governs mostly 
small angle scattering while both Raman type interac- 
tions give rise to mostly large angle scattering. Mathe- 
matically this is due to the fact that diagonal elements 
of the Franck-Condon factors are sharply peaked around 
axis defined by pump wave vector fco, while off-diagonal 
elements favor off-axial scattering for traps of reasonable 
size ~ a few times the resonant wave length Aq. The in- 
teraction strengths of Rayleigh, Stokes, and anti-Stokes 
processes all depend on the spatial distribution of atoms, 
the initial ground state population A'o, and the laser in- 
duced effective dipole interaction strength g. In general 
these different factors compete and complicated pictures 
of light scattering emerge. But for the current model sys- 
tem, we find that the major role is played by the geom- 
etry factor of the system through r/„^m, of which several 
off-diagonal elements are displayed in Fig. ^. 



III. SMALL ANGLE SCATTERING AND 
OSCILLATORY SUPERRADIANCE 



Since Hji is the leading order interaction term, we will 
first consider its effect by discussing the process of small 
angle scattering. In this case, the system dynamics is 
determined essentially by _ff = Hq + Hr. Using the 
property 77„,„(/-i:) = r]n,n{~K), and taking t) = -^7^ 
with 7j real, we obtain the following Heisenberg operator 
equations of motion 



dt 



d 
dt 



dk"fj^'qnniK)Pj^Cn, 



(5) 



where we have defined auxiliary quadrature operators 
P^ = — bt, = ^fc + ^^^d the form function oper- 
ator F{K) = J2n iln,n{K)Nn. Cn is a slowly variant form 
of Cn- These equations and the auxiliary operators are 
similar to those in the theory of collective atomic recoil 
laser (CARL) Since Pt = -Pj:, we see that the 

atom number operator cJjCn is time independent, a testi- 
mony of motional state number conservation in Rayleigh 
scattering. The system of Eq. is integrable and its 
solutions allow for the calculation of the number of scat- 
tered photons 



sin a;fc</2 
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We note that the amphtude of scattering hght inten- 
sity depends on atomic number fluctuations, the variance 
a'^{Nn) = (iV^) — {Nn)^. Depending on the gas tempera- 
tures, this fluctuating part could dominate over the first 
coherent part at larger scattering angles. At very low 
temperatures, assuming a macroscopic condensed atomic 
population in the lowest motional state n — (0,0,0), we 
find 

{bib,) = 27S (^^^^) ' VoAKfiN'o + mo + eN^)], (7) 

where ^ = is for a classical gas, and e = (1) stands for 
a coherent (Fock) state of the Bose-Einstein condensate 
with an average of Nq atoms. The prefactor in the above 
two equations involves siniujkt), a term from the time 
dependent spectra of the square pump pulse. 

At higher temperatures, it was previously calculated 
that 1], 

n 

n ^ ^ 

where /3 = l/fc^T is proportional to the inverse tempera- 
ture and the trap ground state size is Ua — \J {l/2MuJa) 
in direction a = x,y, z. We can readily see that the 
incoherent part contributes mostly at higher tempera- 
tures while the coherent part is more effective at lower 
temperatures. In order to suppress incoherent large an- 
gle scattering we need to satisfy ^K^a'^f^LUa ^ 1, «-e- 
EaP/2 1, where Ea = K^/2M is the one dimensional 
recoil energy. We emphasize that this result is indepen- 
dent of the shape of trapped atomic sample. In terms of 
recoil temperature Tr — /2MkB, a sufficient condi- 
tion is T <C 0.5Tr. At such low temperatures coherent 
Rayleigh scattering dominates. In order to suppress co- 
herent scattering, we can now take advantage of trap size 
parameters. For a cylindrical sample with L (W) the long 
(short) axis length as in Fig. ^ the coherent scattering 
is controlled by 

l%.0(fc - fco)l' = 6-'='^' -n^ 9cos^ ^^-k-W- sin^ 9«in= ^ ^-W^ (k cos O-ko)- ^ (^q) 

where z axis is chosen to be along the pump field di- 
rection and 9, (p are polar and azmithual angles of the 
scattering direction. Typical geometries of current traps 
are L = W ^ \Q (/im) for spherical traps and L — 
\{)W ^ 200 (/im) for cigar shaped traps. Putting these 
into Eq. (M), we see that the cigar shaped geometry 
with L ^ W is more effective in suppressing the overall 
coherent Rayleigh scattering. 

Concluding this section we note that under optimal 
conditions, it is possible to significantly suppress the co- 
herent Rayleigh scattering. Therefore effects of the other 
wise higher order Raman processes can be made domi- 
nant. As reasoned before we can also ignore contributions 
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from anti-Stokes terms at sufficient low temperatures and 
for short pulse excitation. Thus we shall first develop a 
simple model considering Hs as the only dominant mech- 
anism in describing the directional and exponential su- 
perradiance. A more complete treatment including anti- 
Stokes processes will be considered in Sec. V using a 
generalized atomic side mode formalism. 



IV. LARGE ANGLE SCATTERING AND 
EXPONENTIAL GAIN BEHAVIOR 

In this section we consider large angle Raman scat- 
tering processes when small angle Rayleigh scattering is 
suppressed. In the low temperature limit when all atoms 
are condensed into the ground motional state, we can ne- 
glect the ground state atom number fluctuations and ap- 
proximate Co — \/Nq as a classical field. The subsequent 
system dynamics for light scattering can be described by 
the Hamiltonian 

= ^ Enc\cn + / dkujkbtbf: 

n>0 '' 

g*ik,t)Y,Vn,o{k-ko)cibi + h.c. . (11) 

n>0 

In this limit, further insight into this problem can be 
obtained by introducing atomic side mode operators 

n>0 

Physically these are wave packet operators in the ground 
state manifold due to absorption of a pump photon fol- 
lowed by emitting a spontaneous photon, resulting in a 
net momentum transfer of q. Because momentum con- 
servation is in general violated during photon absorption 
and emission among two distinct initial and final mo- 
tional states. These wavepacket operators are the best 
compromise one can construct to reflect the momentum 
conservation law for starting in the ground motional state 
before the absorption and emission cycle. Mathemati- 
cally, they satisfy the following commutation relation, 

[f{q), Pip)] = vooiq-p) ~ voo{q)voo{~p) = 'Dqip}- (13) 

Of particular interest is the special case of [f{q), fHo)] — 
1 — TjQQ^q). These operators can thus be visualized as a 
deformation on the Weyl-Heisenberg algebra [11(4)] of the 
original c„ operators. On the other hand, if we were to 
keep Co as an operator, we would need to consider Xno — 
cJjCo operators We find that c„ can be realized as a 
bosonic representation of transition operators Xno , which 
obeys a U(3) Lie algebra. A similar deformation on the 
U(3) could also be proposed. More generally, a class of 
side mode operators from an arbitrary motional state, 
not limited to n = could also be considered as in Sec. 
V. They represent collective recoil modes corresponding 
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to various transitions from any motional state to higher 
or lower motional states with a net recoil momentum q 
during the photon absorption and emission cycle. 

Neglecting the Doppler effect but keeping the recoil 
energy term ujii{K) = K'^ /2M , we can show that 

EnmA<l)Cn ~ [Ea + LOB{q)]f{q)- (14) 

n>0 

The equations of motion for the operators are then given 
by 

4/(9) = ^R{q)f{q) + / dkg*{k,t)b\V^{k^ - k), (15) 



dt 

d 
'Jt 



i-Jk^^kbj:+^/N^9*(k,t)fHko-k). (16) 



We can proceed with standard technique to eliminate the 
scattered field modes by substituting the formal solution 



ofEq. (16) into Eq. (hj). This yields 



-fiq) + iLUR{q)fiq) = Fl+No J dk g* {k,t)V^{k^ ~ k) dr f {k^ ~ k,t - T)g{k,t ~ T)e'^-\ (17) 

where the Langevin noise operator F^lt), representing 
the effect of vacuum fluctuations through the initial scat- 
tered field operators &j:(0), is introduced as in 

FUt) = J dkg{k, t)Vg-{ko ~ k)b^{Q)e~''^''K (18) 

This noise term is responsible for triggering the super- 
radiant emission from the gas sample. It is also needed to 
satisfy thermal dynamic requirement of fluctuation and 
dissipation theorem. Since its magnitude is of a lower 
order in ^/Nq, it can be neglected in a semi-classical de- 
scription of the growth behavior of a small signal gain as 
in Ref. In the Markov approximation and with typ- 
ically small recoils [ujn^q) 0], we can ignore the slow 
time dependence of g{k,t) to obtain 

d f 

-f{q) = nNo J dk\g{k,0)\^Vi^iko - k)f{ko ~ k)5{ojk)- (19) 

This expression can be further simplifled by noting that 
Dqiko — k) limits the scattering to directions around 
the end firing modes kg — q as illustrated in Fig. |^. 
We then obtain a simple exponential gain behavior 
for atomic side mode operator according to f{q,t) = 
exp {NoG^t/2)f{q, 0) with the gain parameter given by 

G,-=ttJ dk\gik,0)\^Vg-iko - k)5{iUk). (20) 

It is interesting to note that this result is essentially the 
same as obtained in Ref. ||] except for the difference in 
2?q-(p). In fact, the condensate shape function was defined 
in Ref. |^ as Pq{k) = J dr |0o(^)P exp[— i(fc — kQ + q) -r]. 
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It is related our function defined here by 'D^{kQ ~ k) = 

Pq{k) - Pij{ka)pt){k). 

A more rigorous treatment would include explicitly 
the role of atomic quantum statistics. Treating A^o ^ 
an operator, and defining the number operators of re- 
coiled atoms as nq{t) = {f^{q,t)f{q,t)'), we find n^{t) = 
Gq{t)nq'{0), where the growth function is now defined to 
be Gq{t) = (exp (G^iiVo)). In normally ordered form it 
can be written as 

oc 

g,^{t) ^ (e*«f^°) = ^[e*«'^ - l]-'"(4Vo)/j!. (21) 

Clearly, for different quantum statistical distributions, 
this implies different growth behavior. If we choose the 
initial condensate to be a coherent state in the motional 
ground state with an amplitude a such that (cJcq) = 
|ap — (Nq), we obtain 

g,^(t)=exp[(eG?*-l)(7Vo)]. (22) 

While for a Fock state distribution, the same result as 
in the semi-classical case applies. We note that at early 
times {t <C both growth curves give identical re- 

sults irrespective of the atom number statistics of the 
condensate. At longer times, however, an initial conden- 
sate in a coherent state causes side modes to grow faster 
than an initial Fock state condensate. In Fig. |^, we com- 
pare the different growth curves for both cases. Given 
the same number of condensate atoms, the corresponding 
super-radiant pulse is then shorter for a coherent state. 
For the recent MIT experiment Q , it was estimated that 
with a laser intensity of 10 — 100 (mW/cm^), < 10^*. 
Thus for all experimental observed duration, the growth 
curve is the same irrespective of atom number statistics. 
We note that if the MIT experiment were operated with 
a higher pump laser intensity, condensate atom number 
fluctuations could be probed. In contrast to the small 
angle Rayleigh scattering where number fluctuations ap- 
pear as amplitude fluctuations, the large angle Raman 
scattering studied in this section carries information di- 
rectly related to condensate number fluctuations. 

V. SEQUENTIAL SUPERRADIANCE 

In the previous section we considered short pulse 
Stokes Raman scattering in which processes starting from 
the motional ground state n = (0,0,0) dominate. As 
a result, the momentum distribution of atoms, sharply 
peaked around the center-mode p = initially, is mod- 
ified by the appearance of side mode peaks around p = 
ko±ke, where ±fee denote the two end firing modes, i.e. 
emissions along the two ends of an elongated condensate 
as in Fig. |^. This situation is reminiscent of earlier stud- 
ies of superfluorescence from an extended and inverted 
medium [p6|-p2|. When the Fresnel number of the sys- 
tem, defined as — tt{W/2)^ / XqL, is of order unity, a 
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description of the emission can be made in terms of the 



two end firing modes |16|. For the MIT experiment with 
W ^20 (fim), L - 200 (/im), and Aq 0.6 (/xm), the con- 
dition 1 is indeed satisfied. With the pump incident 
along the narrower side of the condensate, possibihties 
exist for mode couphngs into similar end firing modes. 
This in turn causes recoiling atoms to couple with side 
modes of p — ko ± k^, and even higher order side modes 
if the pump stays on for a long period of time. Since 
flo,m{k{) ± ke) is peaked around certain motional state 
mi, this will be reflected in the measured density profile 
as effective couplings populate state mi from condensate 
atoms in state m = (0,0,0). With the wave packet for- 
mulation, a breadth of ground motional states centered 
around mi is affected, and it is termed collectively as 
a side mode to the original condensate. This physical 
picture is further illustrated by the side mode lattice as 
given in Fig. where the most important coupling terms 
from the Hamiltonian Eq. (^) have been selected. We 
can then truncate the number of nodes involved in the 
lattice depending on the duration of the pulse to study 
light scattering dynamics. We introduce generalized side 
mode operators as follows, 

fn{q) = ^ Vn,7n{Q~ P)fva{p), (23) 
m 

with fn{0) — Cn- Their commutator is evaluated to be 

We can then transform the system Hamiltonian into these 
side modes to obtain, 

H = J2[En+UJR{q)]fl{Q)fniq} 

n 

+ X!^^ ■ ^n,™/«(9)/™(9) + / dkujkbibj: 

n.m 

+ dk[g*{k,t)fl{q-k + ko)blfn{q) + h.c.]. (25) 

In the second term, vr — q/M is the recoil velocity and 
Pn,m is the CM momentum matrix elements between 
the Fock basis |n) and |m). This term couples /„ and 
fn+i and in effect describes the hopping between near- 
est neighbor motional states. In our discussion below we 
will neglect this nearest neighbor motional coupling as it 
is_ small in the short time scale of the MIT experiment 
10[. Keeping only the two end firing photon modes 



around ±fce, we can use Eq. (23) to express side modes 
around them according to 

fn{q + fco - fc) = ^ r]n,miTke " k)fm{q + ^0 ± ^e)- (26) 



This implies that by simply examining the dynamics of 
the two end firing side modes, we can also gain valu- 
able understanding of the behavior of other side modes 
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around them. We further simphfy the problem by tak- 
ing only the diagonal Franck-Condon factors, justified by 
the reasonably small Fresnel number T . This allows us 
to use /„((f + fco - ^) ~ Vn.niTke " k) fn{q + ko ± he) as 
(=Ffce - k) 6nm for modes near the two end firing 
ones with k « T^e- Since only /o(0) is initially occupied, 
this Hamiltonian then couples side modes with fo{q) for 
q — 0,kQ± fee, 2fco, 2(fco ± fee),' ■ ■ etc. through an infinite 
hierarchy of equations of nearest neighbor coupling on 
the triangular lattice as in Fig. ||. 

We now discuss effects of the second order side modes. 
Since the central side mode at 2fco is coupled to two first 
order side modes at ko±ke, it will grow faster than other 
second order ones at 2(fco ± ke). Therefore, as indicated 
in Fig. 1^, we close the system of coupled nodes by con- 
sidering the 4 lattice nodes at 0, ko zL k^, 2fco, which are 
connected with solid lines. Pulses with longer duration, 
however, would result in populations grow at higher or- 
der lattice nodes. After free expansion on turning off the 
trapping potential, this particular lattice structure is in 
fact directly observed in the MIT experiment Q. The 
effective Hamiltonian is now given by 

H=J2^Ri'l)foi<l)fo{q)+ I dkL^kbib^ 

Q 

+ J2 [9*At)Blfl{q~ek, + k^)h{q) + h.c], (27) 

where q runs over lattice sites 0, fco ± k^, such that we 
have a truncated problem on the first diamond. The lat- 
tice only runs in one direction with a positive fco because 
of the plane wave excitation from one side. The first term 
of Eq. ( p7[ ) is due to the recoil shift, and can be elimi- 
nated by transforming to an interaction picture. g{k,t) 
is a slowly varying function of k and is replaced with its 
value at k — ke and taken out of the integration over k. 
The emission photon wave packet operator is defined as 
Bf: = J dkr]ofi{k — eke)hj^ for k near ±fce. The evolution 
of its corresponding intensity naturally gives the photon 
scattering distribution averaged over repeated single-shot 
experiments. As a collective field operator, takes into 
account the multi mode but directed (end firing) nature 
of the scattered field [Q. Although mathematically one 
obtains 

(28) 

We take [Be, B\,\ = S^^i since the choice of keeping only 
end firing modes constrains k to be around zt/cg- Thus 
one simply has B^ ~ . Other modes around the end 
firing ones in B^ only contributes to a renormalization of 
the coupling constant g which we take to be phenomeno- 
logical. Terms involving /[2(fco ± /ce)] are ignored in this 
study based on arguments of short pump pulses and low 
atom number populations. 
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It is now useful to introduce a more concise notation 
ao, a±, and 02 for /o(0), fo{ko±ke), and fo{2ko). We also 
treat these operators as commuting with each other as an 
approximation to Eq. (^) in the limit of large |g — p|. 
Their Heisenberg operator equations can be derived. It 
turns out their dynamics is more transparently expressed 
in terms of the population operators I± = i?±S±, = 
alue, and coherence operators R^^' = ala^'. We find 



''l 


= igRo^B± -+ 


igR±2B± + h.c, 




= -igRa-B+ 


— igRo+B^ + h.c. 




= igRo±B^ - 


igR±2B± + h.c, 




= igR+2B+ + 


igR^2B- + h.c. . 



(29) 



In deriving this and other equations to follow, we have 
consistently used an operator ordering with atomic oper- 
ators always to the left of all photon operators. A careful 
analysis Eq. (|2^) reveals that the following two conser- 
vation laws are observed 

7Vo(t) + N^{t) + N+{t) + N2{t) = Ci, 
TVo (t) - (t) + /+ (t ) + /- (t ) = C2 , (30) 

with constants Cj determined by initial conditions. In 
fact, the second conservation law immediately implies 
possibility of sequential superradiance. In the early 
stages of the applied pump when condensate depletion 
is small, scattered light intensity remains low, although 
gradually increasing. Eventually, the rapid decay of the 
condensate populations (TVq N±) sets in and the total 
light intensity starts to increase sharply. The scatter- 
ing losses and absorption then causes the light intensity 
to decay and finally vanish (when A'o empties) while N2 
remains small. The dynamics upto this point is indeed 
equivalent to a system without the presence of N2 term, 
and is simply a parametric amplification process. On the 
other hand, for long pulses with sufhcicnt intensity, the 
now populated N± nodes start to dynamically populate 
node iV2. Thus allowing for an revival of the scattered 
light intensity. 

In the following, we shall be most interested only in 
the population dynamics. Instead of solving the full set 
of Eq. (^9|), we assume equal population distribution 
among the symmetric nodes of Fig. 5|, i. e. treating nodes 
of k[j ± /cg as equivalent. We can then define I = 1+ + 1- 
and Ni — + iV_ . This allows for the consideration of 
an effective set of equations 

= igRoiB + igRi2B + h.c, 

at 

■^Nq = -igRaiB + h.c, 
at 

^Ni = igRoiB - igRi2B + h.c, 
at 
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j^N2=igRi2B + h.c., (31) 

with B now denoting either of B±. The same conserva- 
tion laws Nq + Ni + N2 = Ci and iVo - iV2 + / = C2 
apply. We can proceed to eliminate the scattered field 
operator B from the population dynamics Eq. (^ by 
using the standard technique of substituting in the for- 
mal solution for h-^{t). A Markovian version of closed 
equations will be obtained this way later which allows 
for a direct numerical simulation in terms of averaged 
variables Alternatively, we choose to develop a hier- 
archy of equations for various operator moments first. It 
is illuminating to follow both methods and compare their 
results in the end. 

We now introduce operators X = RqiB, Y = R12B, 
Z = R12R10, and W = i?02-S^. This is needed for a more 
rigorous treatment of operator correlations, a procedure 
similar to the random phase approximation p^ ] . A trivial 
first order decorrelation approximation between matter 
and field (RqiB) = (Rqi) (B) would have neglected too 
much correlation. By forming products involving at least 
4 operators from 6t, 6^, a|, and a^, and making corre- 
sponding decorrelation approximations {IX) = (/) {X), 
etc, we aim for a closed set of equations. Although sys- 
tematic, this approximation procedure is not necessarily 
self-consistent as neglected higher order correlated terms 
may be of the same magnitude of the kept moments. 
There is solid evidence that this is a good approxima- 
tion for super-radiant systems (24j. Our aim is to ob- 
tain a truncated set of equations involving only limited 
number of higher order operator moments which are rel- 
evant to the population dynamics. We start by taking 



the averages of Eq. (31), the right hand side then imme- 
diately motivates the introduction of operators X, Y, Z, 
and W. Upon averaging over their dynamic operator 
equations, even higher order moments in general appear. 
We then follow the decorrelation approximation as out- 
lined above, and only keep factorized products of already 
introduced operator moments. More details of this higher 
order decorrelation approximation can be found in previ- 
ous treatments of superfluorescence . Finally we drop 
the sign (.) for averages and replace g by —ij with 7 real 
to obtain 

^I = -kI + -/{X + Y + c.c.), 
at 

^No = -^iX + C.C.), 
at 

j^Ni = -CNi +-riX-Y + c.c), 

±N2^^{Y + c.c.), (32) 
equations for the higher order operator moments 

= -7iX + 7[iVo(l + Ni) + I{No - Ni) + Z* - W], 

^Y = -ji_Y + j[Ni{l + N2) + I{Ni -N2) + Z + Wl 
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7±Z - 2-f[X*{N2 - Ni) + Y{Ni - No)], 
j^W-2j[IiY-X)^NoY~Xil + N2)i (33) 

as well as equations for their complex conjugates. 

This system ol twelve equations can be compared 
with the Maxwell-Bloch equations describing superfluo- 
rescence from a sample of coherently pumped three level 
atoms |2^,^. Similar to current trapped BEC systems, 
earlier superfluorescence studies also assumed cigar (or 
pencil) shaped gas distribution, albeit with much larger 
volumes. In those earlier studies, the pump pulse is typ- 
ically along the long axis direction of the gas sample 
that is typically about several centimeters long Jl7| , pjt . 
Theoretical analysis included both propagation retarda- 
tion and transverse diffractive effects p5| , p6[ . In the re- 
cent MIT experiment, the sample size is much smaller 
and the pump pulse is along the transverse direction of 
the long axis y. Nevertheless, at higher pump power 
super-radiant pulse shapes from the BEC display multi- 
ple pulses or ringing effects similar to earlier hot gases 
experiments p^ . In this respect, sequential superradi- 
ance may arguably be considered as a temporal analog of 
spatial effects observed earlier although the mechanism is 
clearly different. More recent experiments observed tem- 
poral ringing as an intrinsic property in hot gas super- 
fluorescence . Here in the BEC superradiance system 
ringing can also be understood in terms of the cascading 
structure on the lattice (Fig. ||) as opposed to among 
different electronic levels [ p9[ . 

In the above Eqs. (^) and (|3^), we have introduced 
phenomenological parameters j±, C, and k which are re- 
spectively atomic side mode dephasing rate, decay rate 
of first atomic side modes due to coupling with excluded 
nodes at q = 2(fco ± kg), and the linear loss of scattered 
field in the Maxwell-Bloch equation. Within the time 
scale of interest, losses in N2 due to its coupling to third 
order side modes are negligible. The motional ground 
state-condensate node, is coupled only to the two first 
order side modes thus no dissipative terms appear in the 
equation for iVo. In view of the one to one correspondence 
between the number of atoms in the side modes and their 
corresponding number of scattered photons, we will fur- 
ther assume k — C. In Ref. [|j it was estimated that 
the decoherence time (the decay time of matter wave 
interference) was found to be '--^ 32 — 35 (ps). In our 
numerical simulations, we will thus take dephasing rate 
varies in the range 7^ ^ {2tt) 0.8 — 1.6 x 10* (Hz), while 
C ^ (27r)0.5 X 10"* (Hz). The coherent coupling parameter 
g depends on the pump laser power. Ref. Q reports a 
typical Rayleigh scattering rate ~ {2n) 7 — 700 (Hz) at 
pump intensities ^ 1 — 100 (mW/cm^). Wc thus take 
7/(27r) =5-15 (Hz). The initial condensate number is 
chosen to be A^o = 4 x 10^ [Q for all numerical simula- 
tions. 

In Figures |6|-§, -f± — {2n) 1.6 x 10** (Hz) is used, while 
7i = (27r) 8 X 10^ (Hz) is used in Figs. |l^]l|. The deple- 



dt 
dt 
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tion of condensate atoms is shown in Figs. |6| for several 
different choices of 7/(27r) =5.1, 6.7, 8.3, and 10.7 (Hz) 
with larger coupling rates correspond to faster decaying 
curves. The effect of on the condensate decay can be 
obtained from a comparison of Fig. ^ with Fig. in 
which only two curves with ^/{2ti) — 5.1 and 10.7 (Hz) 
are displayed. Looking at these two curves, we can iden- 
tify four separate dynamical regimes. First, there is the 
linear regime where the condensate can be considered as 
undepleted. Second, there is the super-radiant regime 
where a fast decay of condensate atom number occurs. 
The third regime is a transient one where oscillatory be- 
havior is seen. Finally we have the fourth regime where 
saturation sets in and the condensate decays slowly. Figs. 
^ and ^ show that oscillations are more prominent at 
smaller decoherence rate and higher pump laser power. 
In this case, the linear regime is shortened and super- 
radiant decay becomes faster. 

Figures and ^ display matter wave amplifications for 
the first and second atomic side modes respectively for 
the same parameter set as used in Fig. ^ They should 
be compared with Figs, pi] and ^ where the same pa- 
rameter set as in Fig. |l^ are used . We note that shorter 
and more intense super-radiant pulses are obtained at 
higher pump laser powers. Figs. || and [l3| represent the 
temporal evolution of scattered light intensities for the 
two sets of parameters used for Figs. || and ^ respec- 
tively. They are seen to increase sharply but to decay 
more slowly. At higher laser powers double peak (shoul- 
der) structures are seen while at lower laser powers the 
pulse shape becomes more asymmetric, broad and has 
single peak. With lower decoherence rate the double- 
peak structure is more prominent, and the pulse becomes 
more intense and shorter. Additional peaks in the rings 
are also observed in tail regions. Physically we assign 
the double-peak structure as due to sequential (cascad- 
ing) super-radiant scattering. 

Overall we find that the high order decorrelation ap- 
proximation produces simulation results capable of cap- 
turing detailed dynamics of the super-radiant scattering 
from trapped condensates Apart from the unavoid- 
able choice of introducing and adjusting phenomenologi- 
cal coupling constants and various decay and decoherence 
rates, all reported experimental observations can be in- 
terpreted based on our model |Q. We also find that Z 
and W terms are almost negligible in affecting system dy- 
namics, presumably because they are already of a higher 
order as compared with X and Y . Further research into 
this point is needed for a complete characterization of 
quantum states of scattered photons. 

For completeness, we finally discuss solutions to Eq. 
( ^ ) using the Markov approximation. We formally inte- 
grate the Heisenberg operator equation for &g(t) in terms 
its initial condition &g(0) and a radiation reaction term 
related to emitted field from atoms. The formal solu- 
tion is then substituted into atomic operator equations 
and the standard Markov approximation made such that 
the radiation reaction on atoms become instantaneous. 
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The resulting Langevin equations for atomic operators 
now contain quantum noise terms Q due to bj:{0) and 

&t(0). Averaging over this quantum noise reservoir and 
taking care of operator ordering, we obtain the following 
equations for averaged quantities (again neglecting (.)) 

d 

-^No = -jmNo{1 + Ni), 

j^Ni = 7M[iVo(l + TVi) - N^il + N2)] ~ CNu 

j^No^JmNi{1 + N2), (34) 

with 7m the Markovian coupling constant oc 7^. It also 
depends on reservoir noise spectra width as well as a 
shape factor of the condensate. The above procedure is 
similar to what is used in obtaining Eqs. (|l5|), ([l6|), (17) 



(p^), and (|T9|). 7^/ can be estimated from the Rayleigh 
scattering rate R to be 7m — 3Rft / 8tt (where typically in 
the MIT experimental set up [g , 2 x 10"'') . We again 
will use A'o = 4 X 10^ and the new phenomenological loss 
rate [due to scattering into nodes at q — 2(fco ± fee) in 
Fig. 1) is chosen to be £ = (27r) 4.71 x 10^ (Hz). 

The results for condensate and side mode populations 
from Markovian dynamics Eq. (jsj) are presented in Figs. 
|l4-16. The four sets of curves are gain respectively for 
7M7r27r) = 3.18 x lO'^, 6.37 x lO'^, 3.18 x 10"^, and 
1.2 X 10"^ (Hz), with larger 7m values correspond to 
faster decaying condensates. We see that the transient 
regime with oscillations is visibly absent, while previ- 
ous results from the decorrelation approximation [Eqs. 
( ^ ) and (|3^)] correctly captures this essential experi- 
mental feature. In addition, we find the scattered in- 
tensity displays a two pulse structure in the Markov ap- 
proximation, rather than the ringing shoulder structure 
discussed earlier. We thus come to the conclusion that 
a simple Markov approximation is incapable of describ- 
ing dynamic processes of the super-radiant BEC system. 
In the Markovian limit, the dynamical behavior of scat- 
tered intensity follows / iVo(l + A^i) -f- A^2(l + ^^2)- By 
incorporating quantum fluctuation and choosing appro- 
priate initial conditions, it may be possible to find av- 
eraged intensity distribution from repeated simulations. 
The study in Ref. [^ presented results from a collection 
of single-shot simulations. 



VI. CONCLUSIONS 

In summary, we have presented a thorough investiga- 
tion of the super-radiant light scattering from trapped 
Bose condensates. We have shown that depending on gas 
sample conditions: e.g. density, geometrical shapes, tem- 
peratures, and pump field characteristics, the scattered 
light can exhibit either oscillatory or exponential gain 
behaviors. We have presented a new atomic side mode 
formulation that allows for a useful and simple framework 
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to analyze problems related to scattering from trapped 
particles in momentum space. A cascading structure in 
the side mode lattice is presented that allows for identi- 
fication of sequential super-radiant pulse dynamics. Our 
investigations also point to eminent correlations between 
photons from the two end firing modes similar to the two 
cascading photons from a single atom [ ^0| . Interestingly, 
atomic side modes connected through multiple scattering 
are highly correlated because of their non-commuting al- 
gebra. Even though we simplified that aspect of them 
in accordance with the current experimental designs, fu- 
ture experiments might as well take this as an advan- 
tage of generating strongly correlated high flux of pho- 
tons. We are currently investigating prospects of a super- 
radiant source for correlated or entangled photons from 
our model system. Our study also clarifies similarities 
and differences of the recent MIT super-radiant exper- 
iment from trapped atomic BEG with earlier studies of 
hot gas superfluorescence. We have compared theoretical 
approaches based a Markov approximation with a non- 
Markovian description of the experimental light scatter- 
ing observation. We find that the occurrence of mul- 
tiple peaks in super-radiant pulse reflects the generally 
non-Markovian nature of our system. Our studies also 
sheds light on the role of dephasing in a coherent quan- 
tum system Finally we note that a simple modi- 
fication in atomic side mode definition allows for stud- 
ies of super-radiant emission from a quantum degenerate 
trapped fermi gas. 
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FIG. 1. The geometry of the scattering arrangement from 
trapped atoms. The filled ellipsoid represents the atomic 
cloud with dimensions L and W . Incident pump field comes 
in along the short axis of the trap, while the geometry favors 
the emission to be along the long axis of the trap. Momen- 
tum conservation of the absorption and the subsequent emis- 
sion of photons results in the recoil of trapped atoms into a 
wavepacket state parametrized by q. 
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FIG. 2. The diagram for Rayleigh, Raman Stokes and 
anti-Stokes scattering among the motional states of trapped 
atoms. The solid lines denote pump photons while dotted 
curve lines are for scattered photons. Solid dots denote the 
presence of an atom, and hollow dots denote the absence of 
atoms due to scattering out of certain motional states. 
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FIG. 3. Typical values of the Franck-Condon factor 
\Vo,n{0 \ for C = 0.3, 1, and 3. 
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FIG. 4. Comparison of different growth curves Gqit) for an 
initial motional ground state condensate with different statis- 
tics. Dashed line is for the Fock state and solid line is for 
the coherent state. Semi-classical treatment also predicts a 
growth same with the Fock state result 
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q = 

FIG. 5. The lattice of strongly coupled atomic side modes. 
Solid lines show the diamond shaped unit cell with cascade 
dynamics. 
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FIG. 6. The decay of condensate population. 
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FIG. 7. The atomic population dynamics for the first side 
mode(s). 
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FIG. 10. The decay of condensate population. 
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FIG. 11. The atomic population dynamics of the first side 
mode(s). 
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FIG. 12. The same as in Fig. 8 but for the second side 
mode. 
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FIG. 13. The intensity of scattered field intensity around 
one of the end firing modes. 
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FIG. 15. The evolution of the first side mode(s) within a 
Markov approximation for scattered light. 
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FIG. 16. The growth of second side mode within a Markov 
approximation for scattered light. Curves are magnified for 
better view. 
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